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Abstract. We introduce a meta-population version of models of asymmetric 
exclusion models, consisting of a spatial arrangement of patches. Patches are of 
a specific size, indicating the maximal number of particles they can hold. We 
use an expansion in the inverse patch size to calculate the spectral properties 
of fluctuations in such systems. This provides a systematic derivation from first 
principles of effective Langevin descriptions discussed in the literature. We apply 
our approach to the totally asymmetric simple exclusion process, to variants 
with an overall constraint on the total particle number and to a two-species 
exclusion model. The theory provides semi-analytical results, these are confirmed 
in numerical simulations, and give good approximations to conventional exclusion 
models. These are recovered when the patch size is set to unity. 
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1. Introduction 

The asymmetric exclusion process (ASEP) is one of tlie most studied models in non- 
equilibrium statistical physics. Originally introduced in the context of molecular 
transport [l}|5], it has found wide applications not only to model biological systems, 
but also to the modelling of pedestrian motion and the formation of traffic jams 1 6 lO] . 
Most ASEP models describe a chain of cells, each of which can either be vacant or 
filled by one particle. Particles are injected stochastically at one end of the chain and 
then propagate from one cell to the next according to a stochastic rule. Crucially they 
can only move ahead if the subsequent cell is not occupied, otherwise their motion is 
blocked until the cell ahead becomes vacant. Stochastic ejection occurs at the other 
end of the chain. A comprehensive review of the theory and applications of exclusion 



processes and related models can be found in 11 . In statistical physics exclusion 
processes represent a widely studied class of driven lattice gases 12 - 18 , a variety of 
different methods have been used to characterize their phase behaviour, and to derive 
exact or approximate solutions [8| |THp}fT5l[T7l[l9p2) . 

Formulating a mean-field theory of exclusion processes is relatively straightforward, 
see [2 11 . At the same time it can give accurate insights into the basic phenomena 



displayed by such models. The mathematics required to carry out exact analyses 
on the other hand is intricate, so that it is desirable to develop approaches which 
systematically improve on mean field descriptions, but do not require an overly 
involved mathematical apparatus. Fluctuation effects in ASEP models have for 



example been studied by an effective Langevin description in 23 - 27 . While successful 
in describing the spectral properties of fiuctuations about the mean field theory the 
precise form of these equations is often not derived from first principles. Instead they 
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are formulated phenomenologically, frequently some of the resulting coefficients are 
effectively treated as fitting parameters. 

The purpose of the present work is to discuss a more systematic approach to deriving 
these effective Langevin equations. We start from what we will refer to as a 
'meta-population' version of the totally asymmetric exclusion process. The term 



'meta-population' goes back to Levins 28 who used it to describe 'populations of 



populations'. We use it in the context of the ASEP to describe models in which each 
cell can be occupied by more than one particle. Each cell then becomes a 'patch' 
in which a population of particles can reside. Particles may hop from one patch 
to another according to rules to be specified below. This introduces an interaction 
between the patches, and the aggregate system constitutes a population of interacting 
populations. 

More specifically, we consider a model in which each site can contain up to Q particles, 
where f2 is a fixed positive integer. Instead of being simply occupied or vacant, as 
in the conventional totally asymmetric simple exclusion process (TASEP), sites in 
this extended model are characterized by a filling factor Xi, which can take values 
Xi — 0, 2/r2, . . . , I. Hopping from patch i to i + 1 then occurs with a rate 
proportional to Xi{l — Xi^i), in particular no hopping is possible if the cell ahead is 
fully occupied. Our semi-analytica approach is then based on a systematic expansion 
of the master equation description in powers of fl^^^^. 

While introducing a model with multiple occupancy might seem to constitute a 
significant departure from the original TASEP we believe there is considerable merit 
in studying the meta-population model within a van Kampen expansion picture. This 
is motivated by the following results: (i) In the limit = I the model reduces to 
the standard TASEP; (ii) As we will show, to lowest order in the expansion in n~^/^ 
(i.e. taking the limit ft — )■ oo) the model reproduces the mean field equations usually 
written down for the TASEP, thus the model interpolates between the conventional 
TASEP and its mean field theory; (iii) Taking the expansion to next-to-leading order 
reproduces the Langevin dynamics proposed and studied in 24 ■ 26j ; (iv) Crucially our 
approach is fully controlled, and allows one to state the limitations of the description 
in terms of Gaussian random processes; (v) The relevant coefficients in the Langevin 
(or Ornstein-Uhlenbeck) dynamics are derived from first principles, they do not need 
to be obtained by a fit from simulation data, our analysis may hence also serve as a 
starting point for a better understanding of the 'serious renormalization' of diffusion 



constants and noise strengths reported in 24 . We apply these methods first to the 



basic TASEP and then to what is referred to as the constrained TASEP 25 27 
As a further application we study a two-species exclusion model, to our knowledge no 
calculations based on Langevin approaches have been reported to date for this variant. 

The remainder of this paper is organised as follows: In Sec. [2] we introduce the 
meta-population variant of the single-species TASEP model. A systematic expansion 
in the inverse capacity of the individual cells is the carried out in Sec. [3] Results 
are compared against simulations in Sec. |4j both for the standard totally asymmetric 
simple exclusion process, and for an exclusion process with an overall constraint on the 
total particle number, see also [25)[27]. In Sec. [5] we then address a specific two-species 
exclusion process, before we draw our conclusions and give an outlook on potential 
future research in Sec. [6l 
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2. Model and master equation description 

Our model system consists of L patches (or urns), labelled i = 1, . . . , L. Each urn can 
accommodate up to fi particles, as shown in Fig. [l] We will write ni{t) £ {0, . . . ,n} 
for the number of particles in urn i at time t. The state of the system at any given time 
is fully characterized by the occupation numbers n = (ni, n2, . . . , rti) G {0, . . . , O}^. 
The continuous-time dynamics we will consider is defined by the following transition 
rates (from state n to state n'): 

7o(n'|n) = a{n) ^ '^n;,«i+i, 

Tj(n'|n) = Sn'^,n,~lSn'.^^,ni+i + l, « = 1, . . . , L - 1, 

ri(n'|n) = /3^<5„.,„,.i, (1) 

where dn',n is the Kronecker delta, i.e. Sn',n = 1 for n = n' and = otherwise. 
The first reaction rate, Tq, here describes injections of particles into urn i — 1. Such 
injection is only possible if this urn is not fully filled already (i.e. only if ni < fl), and 
it occurs with a rate proportional to the model parameter a > 0. In order to keep 
the setup sufficiently general to allow for ASEP models with an overall constraint on 
the particle number in the system we assume that a = a(n), where n = X^i'^j- For 
example we will consider models in which a{n) tends to zero when n approaches a 
maximal capacity. The transition rates T!;, i = l,...,L — 1 correspond to moving a 
particle from urn i to i + 1, again this is only possible if the destination urn is not fully 
occupied already. The quantity finally is the rate with which particles are ejected 
from the last urn, i = L, and is proportional to the second main model parameter 
f3>0. 

The time evolution of the probability Pt (n) of finding the system in state n at time t 
is then given by the following master equation 

L 



a,P.(n)=^^ 

2—0 n'^n 



r,(n|n')Pt(n')~r,(n'|n)Pt(n) 



(2) 



Simulations are carried out using the celebrated Gillespie algorithm 29 30], which 
allows one to generate realizations of the stochastic process described by Eq. (|2| . 



3. System-size expansion and calculation of power spectra 
3.1. Deterministic limit 

Deterministic mean field equations can be obtained by multiplying the above master 
equation by n on both sides and then subsequently summing over n. Writing 
Xi{t) = ri^^ {niit)) = X]n"«^*('^) finds 

x.=i;°!i(x)-Tr(x), i=l,...,L, (3) 

where we have written T^(x) — a°°{x){l — xi), as well as T.°°{x.) — Xi{l — Xi^i) for 
i = 1, . . . , L — 1. The term describing ejection from cell L is given by T£°(x) = f3xL. 
We have here introduced a°°{x) for the injection rate in the deterministic limit, where 
X — i^i)- It is important to note that this is a heuristic derivation and 
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Figure 1. An illustration the 'mcta-population' asymmetric exclusion process 
with a capacity of up to f2 = 5 particles per cell. The arrows depict possible 
transitions for the given configuration, hopping from cell 2 to cell 3 is not allowed, 
as cell 3 is fully occupied. 

that an approximation has been made to factorize quantities such as (xiXi^i) into 
{xi) (xj+i). The resuhing mean-field dynamics can then be written as 

±1 = i^i(x) := a(x)(l - a;i) - - X2), 

Xi ^ Fi{:x.) :^ Xi^i{l ~ Xi) - Xi{l - Xi+i), i==2,...,L-l 

XL ^ Fl{x) -.^ XL^i{l- xl)- l3xL. (4) 

Some explanation regarding the nature of these equations is here appropriate. While 
these equations are easy to write down and intuitive in their interpretation it is a- 
priori unclear what exactly the underlying approximations are that have been made 
to arrive at them. The above mentioned factorization for example is a consequence 
of systematically neglecting all fluctuations and e.g. replacing the distribution of 
realized values of by a delta function at its mean value, i.e. Pt(n) — ^ d{n — {n{t))) 
The expansion technique in powers of fi"^/^ we will discuss below can be useful to 
understand these issues. The above deterministic equations are precisely the outcome 
of the lowest order of the expansion, they are therefore exact in the limit — > 00. 
In this paper we are mostly interested in the stochastic dynamics in the stationary 
regime, and as we will see the properties of fluctuations in this regime can be computed 
from the long-time behaviour of the above deterministic equations. Numerically 
integrating these equations, either for constant a, or for a — a{x) in the case of 
constrained the TASEP, we generally find that these equations approach a fixed point 
X* = (xl, . . . ,x*j^) asymptotically. These fixed point values are evaluated numerically, 
and then used for the analysis of stochastic effects, as described below. 

3.2. Leading-order corrections 

Expanding the above master equation in powers of £7^^-^^ allows one to systematically 
characterize fluctuation effects about the above mean-field dynamics. This technique 
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is known as the system-size expansion and goes back to van Kampen 31|| The 
starting point is the decomposition 



= X^{t) + 



(5) 



of the state n of the system into a deterministic part, x(t) and fluctuations ^(t). 
This procedure is straightforward and has been apphed to spatial and non-spatial 



systems for example in 31-40 , the technical details of such calculations are extensively 



described in ^31^ so that we do not report the intermediate steps here. In fact the result 
of the expansion can be written down directly, using the general results decribed in 
the Appendix of 



41 



To leading order in the expansion one obtains the deterministic equations described 
above (Eqs. Q). To first order in Q^^^'^ one finds a linear Fokker-Planck equation of 
the form 



(6) 



describing the evolution of fiuctuations ^ about the deterministic trajectory. We have 
here written di = d/d^i. The matrices J and B are given by 

OF, 



Ji' 



(7) 



and 



B.,, = -T°^A^i,, + [Xr + T°^i]5.., - rrS^+i,,, J - 1, . . . , L, (8) 

and where all other elements of B vanish. We are only interested in the behaviour 
at large times, so all expressions in Eqs. ([t]) and (|8| are to be evaluated at the 
deterministic fixed point. The Fokker-Planck equation Ml) is equivalent to the following 
set of Langevin equations 

j 

where ri(t) is Gaussian white noise of mean zero and with correlations 

{fj,{t)r^,{t'))^B,,6it~t') (10) 
among its components. This can be written as 

Ut) = ~a[x*miit)~[i-x;miit)+^*im2{t) 
+ [i-xi(t)W[x*]Y,m+vi{t) 

i 

Ut) =[i-x*{tm^iit)^^uit)m 

- [1 - x*+,it)Mt) + x*m,+i{t) + 77,(t), i = 2,...,L-l 
Ut) - (1 - xl{t))^L-i{t) - xl_,a{t) - f3Ut) + VL{t), (11) 
(with a'{x) = doi/dx), or in more compact form as 

^ = J^ + »7, (12) 

X We here stress that the expansion parameter in our analysis is the inverse square root of the 
capacity SI of each cell. The resulting theory thus applies in the limit of large, but finite O. The 
number of cells, L, in the system is a separate model parameter, and remains finite throughout. 
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leading to 

[icjl - J] = ^(cj) (13) 

in Fourier space. This can be inverted straightfowardly, and one obtains the power 
spectrum 

(J,{u:)1{lo')) = [{iu:I-I)-'B{-tcol-r)-']^Jiuj + Lo'). (14) 

This set of expressions in principle contains full information about the temporal auto- 
correlations and cross-correlations of the components of ^, and hence describes the 
properties of fluctuations about the deterministic model to the full (in the Gaussian 
approximation we have made truncating the van Kampen expansion after the sub- 
leading term). We will here mostly focus on the fluctuations of the total number 
of particles in the system about the deterministic value x — ^ ^ x* . These are 

given by S,{t) = J^i^ii'^) their power spectrum is obtained as (^^{uj)^{uj')'^ = 
Oneflnds 

(f(c.)^(c.'))=P(c.)<5(c. + c.'), (15) 

where 

Piu) = ^ [(z^I - I)-'Bi^^u;l I^)-'] . (16) 

ij 

We stress that these are the spectra of £,{t), recalling Eq. ([s]) it might be more 
appropriate to refer to them as re-scaled power spectra, given that a factor of f2 has 
been already scaled out in the spectra. As a result the expression in Eq. ( [l6| does not 
depend on 17. In our further analysis this expression is evaluated numericall}J|] and 
will be compared against simulations for a several different variants of the asymmetric 
exclusion process. 

4. Test against simulations 

In this section we will compare results obtained from our analytical calculations against 
numerical simulations of the exclusion process. We first address the standard TASEP 
and then the constrained TASEP. 

4.1. Standard TASEP 

In the standard TASEP there is no constraint on the overall number of particles in 
the system, injection only depends on the occupancy of the first cell, the injection 
rate assumes a constant value a ^ a{n). The dynamics is hence specified by the 
two parameters a and /3, the latter being the constant ejection rate at the end of 
the chain. This system is known to exhibit several dynamic phases, see e.g. [ll] and 
references therein: (i) the so-called high-density phase (HD) at a > /? and /3 < 1/2, 
(ii) the low-density phase (LD) at f3 > a and a < 1/2, (iii) a co-existence phase at 
< a = P < 1/2 and (iv) the so-called maximum current phase (MC) at a,/3 > 1/2. 

§ To this end we have carried out the L X L matrix inversions and multiplications numerically while 
varying u:. This was done mostly for L = 51. At larger chain lengths it might be advisable to 
diagonalize J once and then to work in the corresponding eigenspace. It is then not required to invert 
L X L matrices for all values of oj tested, instead one single matrix inversion is then sufficient. 
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Figure 2. (Colour online) Power spectrum of fluctuations of the total number of 
particles in the system about their deterministic value. Each of the four panels 
corresponds to different combinations of a and /9, corresponding to the LD phase 
(upper left panel), the MC phase (upper right), the coexistence phase (lower left) 
and the HD phase (l ower ri ght). In each panel we show the theoretical prediction, 
obtained from Eqs. l |14|16| , and results from simulations at different cell volumes 
f2, averaged over at least 10 independent runs. The number of cells in the system 
is L = 51. 



We show results for the power spectrum of the fluctuations of the total particle number 
in the four different phases in Fig. [2) In the co-existence and MC phases we find 



algebraic decay of the power spectra, as already observed in 24 , 25 . In the LD 



and HD phases the power spectra again display power-law decay, but modulated by 
damped oscillations, compare again with ,24^,25^ . 

In all four phases we find good general agreement between the analytical predictions 
based on van Kampen's system-size expansion and numerical simulations. Naturally, 
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the theory compares better against simulations when fl is large (see the results for 
= 100 in Fig. [2]), one should here keep in mind that the system-size expansion 
approach is valid in the limit of large, but finite n. Small deviations can never fully 
be eliminated, see e.g. Fig. [3] at small frequencies. We attribute these to either 
fluctuations, limitations in taking the il — oo limit (simulations are necessarily carried 
out at finite fl) and to other numerical effects relating to taking the Fourier transform 
of finite time series. As seen in Fig. [2] the agreement between theory and simulations 
is reasonably good also for = 1, corresponding to the original TASEP, in which each 
cell can hold at most one particle. We re-iterate again that no fitting procedure has 
been carried out at any step of our analysis. Occasionally one finds non-mononotic 
behaviour of the re-scaled spectra as a function of fi, see e.g. the lower left panel 
of Fig. [2] This is presumably due to effects not captured by the expansion to first 
order in (which is inherently based on the assumption that fluctuations are of 

order fi"^/^, and predicts re-scaled spectra hence do not depend on fl). We speculate 
that going to higher orders in the expansion might be able to describe these non- 
monotonicities. We note though that the agreement between theory and simulations 
is less good in the co-existence phase (lower left panel in Fig. [2j the deviation can 
actually amount to up to a factor of two at certain frequencies). This is either due 
to the fact that we have not fully reached the large- limit, or due to the intrinsic 
dynamics of this phase. Mean field solutions typically display a 'kink' in particle 
density (see e.g. (Ts]), and we find that the location of this domain wall can depend 
on initial conditions used to find the fixed point of the deterministic equations, and as 
stated in [l^ on the precise nature in which the limit a — /? limit is taker{|i| We cannot 
exclude that certain phenomena are at work here which we do not fully understand, 
and we would like to limit our conclusions for the co-existence phase to the observation 
that the theory qualitatively reproduces the shape of the power spectrum correctly. 



4.2. Constrained TASEP 



In the constrained TASEP the effective injection rate a depends on the number of 
particles present in the system, and particles enter and leave at fixed rates. Other 



approaches have been considered e.g. in 42 where a gating process was introduced 



We here follow 25 who proposed a model with an overall constraint on the number 



of particles in the system. Specifically we will use 



a{n) — ao tanh 



(17) 



where n is the total number of particles in the system, rim the maximum number of 
particles allowed in the system. The quantity ric is a cross-over parameter, defining the 
precise shape of the sigmoidal function a{n). Adapting this to the meta-population 
exclusion process, in which each cell can hold up to fl particles, the effective injection 
rate is given by 



a(n) — ao tanh 



where Pm is the maximally allowed density of particles, 
appropriately normalized equivalent of n^. 



(18) 

and where pc is an 



II Surprisingly this does not seem to affect the resulting spectrum, which appears to be independent 
of the position of the domain wall. 
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Figure 3. (Colour online) Power spectra of fluctuations of the total particle 
density in the constrained TASEP. Left: /3 = 0.9, oq = 0.3, pm = 0.205, pc = 0.3. 
Right: (9 = 0.3, ao = 0.7, pm = 0.8, pc = 0.7 . In each panel we show the 
theoretical prediction, obtained from Eqs. | |14|16[ l, and results from simulations 
at different cell volumes Q, averaged over at least 10 independent runs. The 
number of cells in the system is L = 51. 



Results are shown in Figjsjfor two sets of the model parameters, chosen to correspond 



to the values used in 25 . As seen in the figure two distinct regimes can again 
be identified, one with an algebraic decay of the power spectrum and another with 
additional oscillatory modulations. The agreement between theory and simulations 
is again good in, with only relatively small deviations at small values of f2. It is 
here appropriate to comment briefly on one difference between our approach and that 
of 



25 . The authors of 25 



phenomenologically derive Langevin equations not too 
dissimilar from the ones we here obtain using a more systematic approach based on the 
cell-size expansion. A further difference between our work and that of 



25 



concerns 

the details of the deterministic fixed point used to compute the power spectra of 
calculations. In 25 conditions are chosen such that a homogeneous density profile 
can be assumed, i.e. x* — x* for all i,j. If this assumption is made the system 
effectively becomes translation invariant, and can be diagonalized in Fourier space 
(the Fourier transform would here be carried with respect to position space, i.e. 
cell numbers). In our calculation we do not make this assumption, but consider a 
general deterministic fixed point x* = (xi, . . . and as a consequence the system 

of L resulting Langevin equations can not easily be simplified. Computing the power 
spectra as detailed in Eqs. ( 14|16 ) therefore requires more computational resources 
than for the approach taken in |25|. Due to this, and to the fact that we simulate 
models of up to f2 = 100 particles per site, results in our figures are mostly limited 
to L = 51, whereas much larger systems are considered in 25 . The different number 



of cells in the system explain the differences between the figures in [25' and ours. 
This applies to the standard TASEP as well as to the model with a constrained total 
number of particles. 
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5. Two-species TASEP 

5.1. Definition 



As a third example we consider a TASEP with two distinct species of particles 18 
For = 1 each cell can then either be empty (0), occupied by a particle of type 1, or 
occupied by a particle of type 2. Particles of type 1 behave the same way as particles 
in the conventional TASEP do, they hop ahead to the subsequent cell if that cell is 
not occupied by a particle of type 1. In other words, particles of type 1 take no notice 
of the presence or absence of type-2 particles. One has the processes 10 — > 01 and 
12 — > 21. Particles of type 2 on the other hand are barred from moving if the cell 
ahead is occupied by a particle of any type, they only move if the cell ahead is vacant. 
Particles of type 2 thus experience an interaction with particles of type 1, one has 
20 — > 02, but no movement in a configuration of the form 21. As one key difference 
to the two models considered before we here choose periodic boundary conditions. 
There is no injection or ejection of particles, the total number of particles of each 
type is conserved throughout. This setup makes the system translation invariant, 
hence simplifying the calculation of the spectral properties of fluctuations. This will 
be detailed below. 

This model is easily generalized to the 'meta-population' case in which each cell can 
hold up to il particles total. If we denote the number of type-1 particles in cell i by 
rii, and the number of type-2 particles in that cell by m^, then we have rii + rui < 
at all times. The meta-population model is the defined by the transition rates 

n,{fl - Ui+i - TUi+i) 



_ m^{Vl - n,+i - m,+i) 

- ^ , (19) 

corresponding to the processes 10 — > 01, 12 — > 21 and 20 — > 02 respectiveljj^] 
Execution of a reaction with rate 7^ 2 will for example lead to the update — 
Ui — l,ni_(_i — rii+i -f l.m^ — > mi + l,TOi+i — TOi+i — 1. In order to capture the 
periodic boundary conditions we wish to address, expressions of the type i + 1 and 
i — I are here to be read 'modulo L\ where L is the total number of sites in the ring. 
A simulation of this process is shown in Fig. |4] 



5.2. Deterministic and stochastic analysis 

Denoting the concentration of particles of type 1 in cell i hy Xi — jf, and that of 
particles of type 2 by = ^ on obtains the following deterministic dynamics in the 
limit n ^ 00: 

it = - Xi{l - x^+i) + (1 - Xi)xi_i, 

i/i = XiUi+i - Xi^iyi + 2/i-i(l - Xi - yi) - yi{l - Xi+i - yi+i). (20) 

Clearly, x* ~ pi, y* = p2 for all i is a fixed point for all pi, p2- This reflects the 
fact that particle numbers are conserved and that the system is translation invariant. 

*i In [18| these different processes are assumed to occur with independent and potentially different 
rates, for simplicity we here focus on the case of equal rate constants. Generalization is 
straightforward. 
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Figure 4. (Colour online) Space-time diagram of the density of particles of type 
2 obtained from a single simulation run of the two-species TASEP. The system 
consists of 512 cells, model parameters are pi = 0.1, p2 = 0.2, each cell can hold 
up to f! = 10 particles. The colour map indicates the density of particles of type 
2, mi{t)/Q in cell i at time t. 



The densities pi and p2 are indeed model parameters of the two-species TASEP, for 

obvious reasons we restrict their choice to non- negative values with pi + p2 < 1. 

It is again straightforward to carry out an expansion in the inverse cell size, calculations 



of this type for spatial systems can for example be found in 40p3 ■ 46 . To first order in 
one again obtains a set of Langevin equations describing the fluctuations about 
the deterministic fixed point. We now have 2L degrees of freedom, Xi,yi (i = 1, . . . , L), 
and we will denote the corresponding fluctuations by <p = (^i, . . . , ^l, Cii • ■ • j Ci)- 
Performing a Fourier transform in both position and time, one finds a linear equation 
of the form 

M(k,uj)ip{k,uj) = rj{k,uj). (21) 

The 2x2 matrix M(fc,a;) and the properties of the white Gaussian noise rj can be 
obtained analytically (details are reported in the Appendix). The power spectra are 
then found from 

Pn(fc,c^)-(ie(fc,^)l'), P22(fc,^) = (|C(fc,^)l')- (22) 

We here remark that the presence of particles of type 2 are irrelevant for the dynamics 
of the type-1 particles. The expression for Pn (w) is therefore exactly the one one would 
obtain for a single-species TASEP on a ring. 



5.3. Test against simulations 

We compare results of the theoretical computations against numerical simulations in 
Fig. [5l As seen in the figure the qualitative agreement between theory and simulation 
is reasonable. We note that the power spectrum in the k — uj plane exhibits a 'rim' on 
which most power is concentrated. In our sign conventions the rim is found at positive 
k and w, and due to the symmetry P{k,ui) = P{—k, —ui) at negative k and to as well. 
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Figure 5. (Colour online) Power spectrum P22{k,u)) of the two-species TASEP. 
Parameters are as in Fig. [4] The left-hand panel shows theoretical predictions 
for P22{k,u)) (obtained from Eqs. | |22| l). The right-hand panel shows results from 
simulations at f2 = 10. Simulation data is averaged over 10 independent runs. 
The colour indicates the magnitude of P22(k,u)) in logarithmic scale (see legend). 



No significant concentration of power is to be expected when fc > and w < or vice 



versa. This is at variance with tlic stochastic travelhng waves observed in 40 , where 
peaks in the power spectrum are found in all four quadrants of the k — w-plane. It 
is here important to stress that the reaction-diffusion equations describing the system 



of 40 are of second order with respect to position, and hence they are invariant 



against reflections x —x. Waves travel in both directions in the model of 40 . In 
the exclusion processes this reflection symmetry no longer applie;|^ particles generally 
travel to the right (the only exception in our model is the reaction 12 — !■ 21, when the 
particle of type 2 effectively hops to the left), as also seen in Fig. |4j 
In order to provide a more quantitative comparison between simulations and theory 
we depict the spectrum P22{k, lo) as a function of uj at several flxed values of k in Fig. 

m 

6. Conclusions 

In summary we have used a meta-population approach to describe several variants 
of the totally asymmetric exclusion process. The term 'meta-population' here refers 
to a setup in which each cell of the underlying spatial structure can be occupied by 
several particles, up to a total capacity of Vt. Carrying out a systematic expansion in 
powers of one derives the deterministic limiting equations in the leading order 

of the expansion, and obtains a set of Langevin equations describing fluctuations 
about this deterministic limit in the sub-leading order of the cell-size expansion. Such 
Langevin equations are not new for the description of exclusion processes, they have 



for example been formulated and used in 24-26 . These existing studies however take 



a mostly phenomenological approach, we feel that the angle taken here provides a more 

+ The terms on the RHS of Eqs. \2Q\ can be written as first-order lattice derivatives of functions 
such as Xi{l - Xi+i) or yi(l - Xi+x - yi+i). 



Fluctuations in meta-population exclusion processes 



13 






Figure 6. (Colour online) Power spectrum P22{k,Lij) of the two-species TASEP 
at fixed values of fc = ^£ with £ = 2, 3, 4 (left to right at the maximum of the 
curves). Parameters are pi = 0.1, P2 = 0.2, the number of cells in the system is 
L = 128. The solid lines show P22{k,u!) as obtained from Eqs. \22\ , markers are 
from simulations (circles: Q = 100, squares: f! = 10, diamonds: f! = 1), averaged 
over 50 — 500 independent runs, depending on the cell size. 



systematic derivation of an effective Langevin dynamics from first principles, and using 
well-controlled expansion techniques. Secondly, the meta-population model allows for 
a smooth interpolation between the mean-field limit and the standard single-occupancy 
ASEP model, ft is also important to stress that results from this expansion do not 
require any fitting parameters, unlike some of the more phenomenological approaches 
considered previously. While we report results mostly for L = 51 cells in the system 
we note that our results are general and apply to any length of the ASEP chain. Our 
theory might therefore also be useful to shed more light on the L-dependence of fitting 



parameters and on the shape of power spectra reported e.g. in 24 25 . We realize 
of course that the meta-population model is a-priori different from the conventional 
exclusion process, in which each cell can be occupied by at most one particle at any 
time. However, as our results show the predictions derived from the sub-leading order 
of the expansion in the cell size agree reasonably well with simulations of the single- 
occupancy modeQ We are therefore hopeful that the approach taken here might 
be useful to investigate other variants of the exclusion processes, for example two- 



lane models (see e.g. 47 ) or models with spatial heterogeneities and/or individual 



slow' sites 26 . Cell-size expansion techniques may also be considered for models 



* It should be noted that this agreement cannot generally be expected to hold very close to phase 
boundaries. 
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with interactions reaching beyond neighbouring cells, such as for example the Nagel- 
Schreckenberg model of vehicular traffic [9]. Work along these lines is in progress. 
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Appendix A. Spectra of the two-species process 

This appendix provides some more details of the calculation of power spectra in the 
two-species exclusion process discussed in Sec. [5] Similar to the standard TASEP 
model and the TASEP with constrained particle numbers one finds a set of Langevin 
equations of the type 

(f^I(f + T] (A.l) 

to sub-leading order of the cell-size expansion. The matrix J is the (2L) x (2L) 
Jacobian of the deterministic dynamics, and rj describes a 2L-component Gaussian 
noise variable. We will first calculate and simplify the Jacobian J, and then address 
the noise. 

Appendix A.l. Calculation of the relevant Jacobian 
To simplify the notation we define fi and gi 

/j = -Xi{l - Xi+i) + (1 - Xi)Xi-i, 

gi = x,yi+i - Xi^iyi + yi_i(l - - y^) - yi{l - x,+i - y^+i) (A. 2) 



as the expressions on the RHS of Eqs. (20). Then one finds 



as well as 



dxi-i 


= 1 


dxi 


OXi+i 




dh 
dyi-i 


= 0, 




dyi+1 


(A.3) 


dgt 
dxi^i 




axi 


dx,+i 




d.n 


= 1 




dh ^ df, 
oyi oy^+i 


(A.4) 



The 2x2 Jacobian in Fourier space is therefore given by 



J(fc) 



(l~Pi)e-*'=-l-Fpie*'= 

-p2e-*'= + P2e'^ (1 - Pi - P2)e-*^- - 1 + (pi + P2)e''' 

cos(fc) - 1 + i(2pi - 1) sin(A;) 

i2p2 sin(fc) cos(A;) — 1-1- i{2{pi + P2) — 1) sin(fc) 

(A 
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Appendix A. 2. Noise correlator 
Writing 

{VaAt)Vb.jit')) = Bab,,jS{t ~ t') (A.6) 

where a, 6 g {1, 2} denotes the two species, and where j, j € {1, . . . , L} stands for cells 
we have 

Bll,^J = (-T* - n)5,,,-i + {T* + T* + T* + Tl)5,, + (-Ti - Ta)^,, 
Bi2,ij = T^Sj^i-i + (— T3* — Tl)5ij + T^Sj^i+i, 
B2i,ij = Tl5j,i^i + (—2^3* — Tl)5ij + r3(5j 

^22,ij = (-r4* - r*),5,- + (T3* + T4* + T* + Tg*)^ + (-r3 - Ts)^,, 

(A.7) 

where we have used the shorthands 

-'l ~ ^iAJ ^2 — ^1-1,1; 

t; = T*,, t; = 7^:3, (A.8) 

and where the asterisk indicates that these rates are to be evaluated at the 
deterministic fixed point. Carrying out a Fourier transform with respect to position 
space (cell number) one finds 

Bn(fc) = (-T2* - Tl)e''' + (Tj- + T* + T3* + T^) + (-Ti - T3)e-^^ 
Bi2{k) = r4*e'" + (-T3* - Tl)5,, + ^3e-^^ 

B2i(fc) = ns.^.-i + (-T3* - t;)% + ^3e-'^ 

B22(/c) = {-Tl - T*)e''' + {T* + n + T5* + T*) + (-Tg - T5)e-*'=. 

(A.9) 

Inserting the fixed point values Tj* = ~ Pi(l ^ Pi ^ P2), ^a* = = P1P2 and 
T5* = Tg* = ^2(1 - /Oi - P2) this becomes 

Sn(fc)=2pi(pi-l)[cos(fc)-l], 

Bi2ik)=2pip2 [cos(fc)-l], 

B2i{k)^2pip2 [cos(fc)-l], 

B22{k) = 2p2{p2 - 1) [cos(fc) - 1] . (A.IO) 

Appendix A. 3. Langevin equation and power spectrum 

The Langevin equation ( |21[ ) is given by 

{iio -I{k))lp{k,uj) = ri{k,uj), (A.ll) 

where we have carried out Fourier transforms both with respect to position and time. 
This can be written as 

M(fc,a;)^(fc,a;) = ^(fc,a;), (A.12) 

i.e. 
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Writing (f = {x,y) and rj ~ {jh,'i]2) the solution of Eq. ( A.12[ ) is given by 

x{uj) \ _ 1 / m22(fc,w) -mi2{k,u) \ ( rji[uj) 



I.e. 

ibj — cos(fc) + 1 — i{2pi — 1) sin(A;) 

— i2p2 sin(fc) ioj — cos(fc) + 1 — i{2{pi + P2) — 1) sin(A;) 

(A 

_ _ _ (A 15) 

y{<jj) J A{k,uj) \ -m2i(fc,w) TOii(fc,a;) J \ ?72(w) ' ' ' 

where 

A = 7711177122 - "T-12TO21- (A. 16) 

From this we have the final result 

(|S(fc,w)P) = 1^ [ITO22P&II + ITO12P622 - 2Re[77l22mi2]&i2] 

(|y(fc,w)|2) = 1^ [r7iiip622 + |m2ip6ii-2Re[77iii7772i]62i] , (A.17) 
where we have suppressed the dependencies on k and uj on the RHS. 
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